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Abstract 

We present a description of the evolution of a polarized Gaussian beam in 
a smoothly inhomogeneous isotropic medium in frame of the eikonal-based 
complex geometrical optics which describes the phase front and the cross 
section of the Gaussian beam using the quadratic expansion of the complex- 
valued eikonal. The linear complex-valued eikonal components are intro- 
duced to describe the influence of the spin-orbit interaction and the defor- 
mation of a polarized Gaussian beam on the propagation firstly in this paper. 
In an inhomogeneous medium, the interaction between the polarization and 
the rotation deformation of the light beam is presented besides the spin-orbit 
interaction, it corresponds to the spin-intrinsic orbital angular momentum 
interaction and makes the correction for the spin Hall effect of a polarized 
Gaussian light beam. 



1. Introduction 

The propagation of a polarized light has been investigated in various 
physical environment [1-16]. In the process of the propagation, there are 
two important observable effects, the spin Hall effect [1-12] and the Berry 
phase [17] which describe the splitting of the propagation trajectory of dif- 
ferent polarized light and the rotation of the polarization ellipse[18, 19,20] 
respectively. The two phenomena are caused by the spin-orbit interaction of 
photons, which describes the interaction between the polarization and the 
extrinsic orbital angular momentum of light [1,4, 10]. 
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The eikonal-based complex geometrical optics[21,22, 23,24] and the so- 
lution of a parabolic-type wave equation have given the description of the 
Gaussian beam propagation in an inhomogeneous medium, and the evolution 
of the Gaussian beam reduces to the solution of the Riccati type ordinary dif- 
ferential equations. One knows that the quadratic expansion complex- valued 
eikonal[22,23] is introduced to describe the phase front and the cross section 
of the Gaussian beam in frame of the eikonal-based complex geometrical op- 
tics. There is the deformation of the Gaussian which is due to the medium 
inhomogeneity[22,23]. we have known there is the spin-orbit interaction for a 
polarized light beam, but whether there is the interaction between the polar- 
ization and the deformation of the light beam? In this paper, this problem 
is investigated using the eikonal-based complex geometrical optics. 

In frame of the eikonal-based complex geometrical optics, One knows 
that there are no the linear components in the complex- valued eikonal[23] 
when the Gaussian beam propagates in the homogeneous isotropic medium, 
but the inhomogeneous isotropic medium can be considered as the weakly 
anisotropic medium for different polarized light in the first-order geometrical 
optics approximation [3]. In order to describe the perturbation in the propa- 
gation, the linear complex-valued eikonal components are introduced in this 
paper. We obtain the result that the propagation of a polarized Gaussian 
beam in an inhomogeneous isotropic medium is affected by the spin-orbit 
interaction and the rotation deformation of the beam. The spin Hall effect 
of a polarized light beam[25] describes the transverse shift of the beam as 
a whole and is determined by the spin-orbit interaction. The rotation de- 
formation of the light beam produces an equivalent intrinsic orbital angular 
momentum[26, 27,28], there must be the interaction between the polarization 
and the rotation deformation of the light beam. This interaction makes the 
correction for the spin Hall effect of a polarized light beam. 

The paper is organized as follows. First, the eikonal equation and ampli- 
tude transport equation which include the first-order geometrical optics ap- 
proximation correction are derived. Then, the linear complex-valued eikonal 
components are introduced to describe the influence of the medium's in- 
homogeneity on the propagation of a polarized beam within the paraxial 
approximation. And the differential equations describing the evolution of 
a polarized Gaussian beam are obtained. Finally, we reexamine the spin 
Hall effect of a polarized Gaussian beam and obtain result that both the 
spin-orbit interaction and the deformation of a polarized Gaussian affect the 
propagation trajectory of the light beam. 
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2. The eikonal equation in first-order geometrical optics approxi- 
mation 



Let's consider the propagation of a monochromatic linearly polarized light 
in a smoothly inhomogeneous isotropic medium with dielectric permittivity 
e(r). There is a small geometrical optics parameter \iqo [21,29], 



A 



Hgo 



«1, 



(1) 



where L ~ \Vefe\~ 1 is the characteristic scale of the medium inhomogeneity, 
A = 2ir/k = 2-nc/u is the wavelength in vacuum, uj is the angular frequency. 
Next, the eikonal equation will be derived in the first-order approximation 
in nao- 

Maxwell equations have the form 



e <9 2 E 

V*(VxE) + -_ 
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then 
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The last term in the left-hand side of Eq.(3) corresponds to the spin-orbit 
interaction of electromagnetic waves [1,10]. In geometrical optics, the prop- 
agation of the electromagnetic wave can be described by the ray trajectory 
in the short-wavelength limit. A ray-accompanying frame (?7i, ?72, t) is in- 
troduced with the unit vectors (e 1 ,e 2 ,l), r is the parameter along the ray, 
which is connected with the ray arc length s by the relation dr = ds/yfe~ c . 
e c = e(r c ), r c = (0,0, r) is the radius vector for the central ray in zero- 
order approximation in \ico- The Lame coefficients of the coordinate frame 
(Vi,V2,t) are [30,31] 



where 



hi — hi — 1, h = h T = \fe 
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The electric field of the electromagnetic wave in the ray coordinates is E = 
Ei&i + E 2 e 2 + = E ± + Enl, En <C E_l. Keeping the terms up to the fi GO , 
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the wave equation for the transverse electric field Ej_ can be written in the 
form 

V 2 E ± + k 2 eE ± + (V lne x V) x E ± = 0, (6) 

the third term in the left-hand side of Eq.(6) corresponds to the spin-orbit 
interaction, it originates from the third term in the left-hand side of Eq.(3). 
One knows a linearly polarized electromagnetic wave is a superposition of 
right-hand and left-hand circularly polarized waves, Ej_ = E + e + + _E_e_, 

e ± = i=(e 1 ±^e 2 ), E± = ±={E lT iE 2 ). (7) 

Because of 

1 x e CT = -iae a , (8) 
we obtain the wave equation from Eq.(6), 

k^ 2 V 2 E a + eE a - ik^al ■ ( V In e x V)E a = 0, (9) 

where a = ±1 denote the wave helicity of right and left circular polariza- 
tions. One can takes E a as wave function of photons and introduces the 
momentum operator of photons — ik$ 1 V, Eq.(9) can be considered as the 
schrodinger-type equation of photons, the Hamiltonian of photons[l,4,10] 
which include the spin-orbit interaction is obtained. Then the equation of 
motion for photons[l-4,6,10] can be derived using the canonical equations, 
which describes the spin Hall effect of photons. 

The eikonal equation which includes the spin-orbit interaction also can 
be derived from the wave equation (9). From Eq.(7), E a can be written in 
the following form 

E a = A(r) exp[ik ilj{r)] exp(-ia^). (10) 

Substituting Eq.(10) into Eq.(9), the eikonal equation and amplitude trans- 
port equation are obtained in the first-order geometrical optics approxima- 
tion, 

(VVO 2 = e + * x V) -1, (11) 

2V^ • VA + AVV = a x VA^) • 1. (12) 
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The second term in the right-hand side of Eq.(ll) and Eq.(12) are the spin- 
orbit correction which corresponds to the first-order approximation in fico- 
The eikonal equation (11) is same as the equation of motion derived from the 
photon's Hamiltonian, and it also describes the spin Hall effect of photons. 
The amplitude transport equation (12) describes the variation of the ampli- 
tude along the ray. If we take A(r) = A(t), the right-hand side of Eq.(12) 
disappear, then V • (A 2 Vip) = 0, this form means the conservation of the 
energy flux along the ray. Next, the propagation of a polarized Gaussian 
beam is investigated based on Eq.(ll) in the paraxial approximation. 

3. The evolution of a polarized Gaussian beam in an inhomoge- 
neous isotropic medium 

In order to describe the propagation of the Gaussian beam in the paraxial 
approximation, two small parameters are introduced, 

Hi = -, H2 = y, (13) 
w L 

where w is the characteristic beam width. Hi and \x 2 describe the angle of the 
beam diffraction widening and the influence of the medium inhomogeneity 
on the diffraction respectively [23]. Three parameters, fj,i , /i 2 and hgo follow 
relation 

Hgo < A*2- (14) 

The propagation of the Gaussian beam in an inhomogeneous medium has 
been investigated with an accuracy of n\ and /j^ [23,30,31], we will keep the 
order in (j,goiA and /igo/A because of the spin-orbit interaction. 
In the paraxial approximation, the eikonal can be written as 

^ = -0 c ( r ) + 5^(771, t} 2 ,t), (15) 

where i/) c (t) is the eikonal on the central ray and 5ip is a small deviation in 
the paraxial approximation. Sip can be written as the following form[23] 

2 2 

<ty(77i,772,T) = ^S 11 (r) + ^-S 22 (r)+r ]l r ]2 S 12 (r)+r ]l S 1 (r)+r ]2 S 2 (r) + S (r), 

(16) 

where Sij (i = 1,2 and S± 2 = S 2 i) are the quadratic complex- valued term 
in rjiTjj. The real and imaginary components of describe the shape of 
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the phase front and the cross section of a polarized Gaussian beam[21-23]. 
The linear complex- valued term, Si, describe the deviation of the center of 
the beam. One knows that the deviation of the center of the beam is due to 
the spin-orbit interaction in the first-order geometrical optics approximation. 
So(t) is a small additional phase in the paraxial approximation. Substituting 
Eq.(15) and Eq.(16) into the eikonal equation (11), 



dip 
dr] 2 



h 2 V dr 



dine dip 

k \ drji dr/2 drj 2 drji 



a ( d In e dip 



(17) 



The dielectric permittivity e can be expanded in a Taylor series, 



e = e c + t] ■ V±e c + -(77 • V±) 2 e c , 

where 77 = rj^ + r] 2 e 2 , V_l = d/dr), V±e c = (V_Le)| r=r 
Eq.(18), the eikonal equation has the form 



(18) 

From Eq.(17) and 



h 2 



dip 
d^ 



+ 



dip 
di] 2 



+ 



dip c dSip 
dr dr 



hh+h2 a f d In e c dip 



k \ drfx dr\ 2 



In zero-order approximation , the ray equation is obtained 

dipc _ 
dr ~ £c ' 



d\ne c dip 
drj 2 drji 
(19) 



(20) 



This equation describes the evolution of the trajectory of the center of the 
beam without the influence of the spin-orbit interaction[29,32]. 
From Eq.(19) and Eq.(20), we have 



h 2 



dip 
dr]i 



+ 



dip \ 



dr] 2 ) 



2 dipcd6ip_ 
dr dr 



d5ip 
1h 



Sr. , — x2 ^ 



2 v»7-V±re c --(»7-V ± e c )^ 



k 



2 a f dine c dip 



d\ne c dip. 



drji dr/2 

Eq.(21) can be expanded in 77, the differential equations for the complex 



dr] 2 dr) 
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parameters and Si are obtained keeping up to the HgoiA an d HGOfAi 

dS 



dr 
dS 



dr 
dS 12 
dr 



L1 +(^i + ^) = «n, (22) 

" 2 + (^12 + SI2) = «22, (23) 



+ 1 S 12 ( 1 S 11 + 5 22 ) = a 12 , (24) 



^ + (*Sn + S 2 S 12 ) - ^ f - °^S U ) = 0, (25) 



, „ „ <7 / dm £ r _ <9 In e 



^ + (SlSl2 + SA2 ) - Wo ^S 22 - ^S U ) = 0, (26) 
a5 °+i(S? + S 2 2 )-^(^S 2 -^ Sl )=0, (27) 



where 



dr 2 V 1 A ' 2k V drn dr) 

_ 1 d 2 e c 3 <9g c <9e c 
JJ 2 drjidrij 4e c drji drjj ' 

Eqs.(22)-(27) describe the evolution of a polarized Gaussian beam in the 
smoothly inhomogeneous isotropic medium, the parameters describe the 
influence of the medium inhomogeneity on the propagation of the beam. 
These equations can be divided into three groups. The first group is Eqs.(22)- 
(24) which describe the evolution of the shape of the phase front and the cross 
section of a polarized Gaussian beam[23]. This group of equations has been 
obtained in [23] and is independent of the other equations. The second group 
is Eq.(25) and Eq.(26) which describe the transverse shift of the center of 
a polarized Gaussian beam. The third group is Eq. (27) which describes the 
evolution of the additional phase due to the deformation of the beam. 

We have known that Eqs.(22)-(24) keep the terms in fij and Eqs.(25)-(26) 
in ncoul- The spin Hall effect of a polarized Gaussian beam[25] describes 
the transverse shift of the beam as a whole, and the transverse shift is ac- 
cord with the spin Hall effect of photons in the geometrical optics. In next 
section, the transverse shift will be reinvestigated, the interaction between 
the polarization and the rotation deformation of a polarized Gaussian beam 
will be presented besides the spin-orbit interaction. 
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4. The correction of the spin Hall effect of a polarized Gaussian 
beam 



The complex-valued eikonal components can be written as S(t) = S r (t) + 
iS^r), S r (t) and S^fY) are the real and imaginary components. Eq.(lO) has 
the following form 



exp[ik ip(r)] = exp[ik ip c ]exp 
= exp[ikQip^\exp 



ik o ( y Su + ~2 S22 + V^Sw + 771^1 + V2S2 



-k (f^! + f^ 2 + V1V2SI2 + msl + msl) 



exp[ik Q i(j c } 



exp 
exp 

where d\ and g?2 are 
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D 22 
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d 



R 



si 

Si?' 



d 



5/ 



(30) 



cif and df are the transverse shift of the center of the wave phase front and 
the cross section[33]. From Eq.(29), one knows that the transverse shift [25] 
don't change the shape of the phase front and the cross section of the beam, 
but the deformation of the phase front and the cross section can affect the 
transverse shift of the light beam. This transverse shift represents the spin 
Hall effect of a polarized light beam[25]. 

By virtue of Eqs. (22-26) and Eq.(30), the evolution equations of the trans- 
verse shift are obtained, 



ddj 



1 dS( S( 
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Because of d[ = S^/S^, Eq.(31) is written as 
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From Eqs.(33)-(36), the transverse shift in the plane (771,772) are 



(40) 



2k J K CJ 2k J\ drn dr]2 
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d ' = a; / (1 x v '^»* + ^ / (i^rst - - ejrsf ) Sl2 + ° ' 

(41) 



where R and O 1 are 



O r = 1 (Of ei + Oie 2 )dr, (42) 
' = y (Of ei + 0?e 2 )dr. (43) 



In Eq.(40) and Eq.(41), d = d 1 e 1 +d 2 e 2 ,d i = d^. From Eq.(40) and Eq.(41), 
the transverse shift of a polarized light beam is different from photons. In 
the right-hand side of Eq.(40) and Eq.(41), the first term corresponds to the 
transverse shift due to the spin-orbit interaction, it describes the influence of 
the interaction between the polarization and the trajectory of the propagation 
of the light beam. The second term is closely related to the quadratic terms 
512. The term S± 2 describes the rotation of the cross section of the beam in 
the plane (771, 772) [23,30,31]. The rotation of the beam produces an additional 
equivalent intrinsic orbital angular momentum[26, 27,28], so the second term 
corresponds to the spin-intrinsic orbital angular momentum interaction [34]. 
The third term describes the influence of the deformation on the transverse 
shift. One has known that the propagation of a polarized light (or photons) is 
influenced by the spin-orbit interaction, but a polarized light beam is different 
from photons because of the spatial energy distribution of the light beam. 
This result means that it is not enough if we only consider the influence of 
the spin-orbit interaction for the propagation of a polarized light beam in an 
inhomogeneous isotropic medium, the rotation deformation of the beam also 
is important. This is the main result in this paper. 

The imaginary components of Sij describe the amplitude(energy) spatial 
distribution of the Gaussian beam. From Eq.(40) and Eq.(41), the splitting 
of the center of gravity for the linear polarized Gaussian beam is obtained 

«• - (d <) Md ')- - i/(.xV.n^ + i / - ^|) SUr. 

(44) 

Here, the first term describes the splitting due to the spin-orbit interaction, 
the second term corresponds to the splitting due to the rotation deformation 
of the Gaussian beam. This is a significant difference between the light beam 
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and photons. 

Although the analytic solutions of the differential equations (22)- (27) can- 
not be obtained, we give the numerical simulation when the linear polarized 
Gaussian beam propagates along the helical ray in the cylindrical symmetry 
medium[23]. Fig.l shows the evolution of the splitting for the center (central 
ray) of the linear polarized Gaussian beam and the propagation trajectory 
of photons. Because of the deformation of a polarized Gaussian beam, the 
relative shift between the right-hand and left-hand circularly polarized Gaus- 
sian beam is slightly different from photons. Although the difference of the 
relative shift is very small, one also should pay attention to the spin-intrinsic 
orbital angular momentum interaction because of the deformation of a po- 
larized light beam. 

5. Conclusion 

The propagation of a polarized Gaussian beam in an inhomogeneous 
isotropic medium has been presented using the description of the eikonal- 
based complex geometrical optics in this paper. The linear complex-valued 
terms are introduced to describe the perturbation of the center of a polarized 
Gaussian beam due to the inhomogeneity of medium and the deformation 
of a polarized light beam. Because the rotation deformation of a polarized 
Gaussian beam produces an additional equivalent intrinsic orbital angular 
momentum, the interaction between the polarization and the rotation de- 
formation of a polarized light beam is presented. This interaction can be 
considered as the spin-intrinsic orbital angular momentum and makes the 
correction for the spin Hall effect of a polarized light beam. 
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Figure 1: Numerical result for the splitting of the linear polarized Gaussian beam and 
photons which both propagate along the helical ray in the cylindrical symmetrical inho- 
mogeneous medium. The red line and the blue line describe the evolution of the splitting 
of the trajectory of the beam and photons respectively. The parameters are as follows: 
the angle between the tangent to the helical ray and the cylinder axis, the dielectric per- 
mittivity e — £o — P 2 /L 2 , £o = 1, L = 200A, the initial beam waist radius wq = 10A, the 
initial phase front is flat r — oo. 
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